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ABSTRACT 

The Meteoroid Input Function (MIF) model has been developed with the purpose of 
understanding the temporal and spatial variability of the meteoroid impact in the atmosphere. 
This model includes the assessment of potential observational biases, namely through the use 
of empirical measurements to characterize the minimum detectable radar cross-section (RCS) 
for the particular High Power Large Aperture (HPLA) radar utilized. This RCS sensitivity 
threshold allows for the characterization of the radar system’s ability to detect particles at a 
given mass and velocity. The MIF has been shown to accurately predict the meteor detection 
rate of several HPLA radar systems, including the Arecibo Observatory (AO) and the Poker 
Flat Incoherent Scatter Radar (PFISR), as well as the seasonal and diurnal variations of the 
meteor flux at various geographic locations. In this paper, the MIF model is used to predict 
several properties of the meteors observed by the Middle and Upper atmosphere (MU) radar, 
including the distributions of meteor areal density, speed, and radiant location. This study 
offers new insight into the accuracy of the MIF, as it addresses the ability of the model to 



predict meteor observations at middle geographic latitudes and for a radar operating 
frequency in the low VHF band. Furthennore, the interferometry capability of the MU radar 
allows for the assessment of the model’s ability to capture infonnation about the fundamental 
input parameters of meteoroid source and speed. This paper demonstrates that the MIF is 
applicable to a wide range of HPLA radar instruments and increases the confidence of using 
the MIF as a global model, and it shows that the model accurately considers the speed and 
sporadic source distributions for the portion of the meteoroid population observable by MU. 

1. INTRODUCTION 

The majority of extraterrestrial particles entering Earth’s atmosphere originate from the 
Sporadic Meteoroid Complex (SMC) and are in the mass range of 10" 11 to 10" 4 g (Ceplecha et 
ah, 1998; Baggaley, 2002; Williams and Murad, 2002). These sub-millimeter sized 
meteoroids enter the atmosphere at speeds between 1 1 and 72 km/s (Close et ah, 2002) and 
subsequently collide with air molecules, causing heating and ablation of the incoming 
particles, which releases both neutral and ionized atoms from the meteoroid as well as ionizes 
the surrounding atmosphere. As a result, meteoroids deposit large amounts of metallic atoms 
(e.g. Si, Na, Fe, K, and Ca) in the Mesosphere and Lower Thennosphere (MLT), mostly over 
the altitude range of 70 to 140 km (Janches et ah, 2003; 2009; Vondrak et ah, 2008), and 
provide the dominant source of metals associated with many mesospheric layer phenomena 
(Plane, 2003). Moreover, similar ionized layers have been observed in other planetary 
atmospheres (e.g. at an altitude of approximately 90 km on Mars, 120 km on Venus, and 550 
km on Titan; Kliore et ah, 2008; Patzold et ah, 2005, 2009; Withers et ah, 2008). 

There are various types of ionization regions associated with the meteoroid, collectively 
referred to as the meteor, which are created throughout the particle’s atmospheric entry 



process (Janches et al., 2003). In particular, there is a relatively small plasma region 
surrounding and moving with the meteoroid, known as the meteor head echo (Janches et ah, 
2000; Close et ah, 2002). High Power Large Aperture (HPLA) radar systems efficiently 
detect this target and have been extensively used over the past decade for various 
investigations that address issues related to plasma physics, space and solar system sciences, 
and astronomy (Janches et al., 2008). These instruments are capable of determining the 
meteoroid’s range and range rate with high precision, and in some cases, the object’s position 
and velocity vectors may be calculated, thus allowing for the estimation of orbital properties 
(Chau and Woodman, 2004; Sparks et al., 2010; Kero et al., 2011). HPLA systems are also 
effective at detecting the evolution of these quantities along the meteor path, and thus, they 
enable the study of various physical processes regarding the meteoroid mass loss mechanisms 
through either the dynamic or scattering behavior of the meteor (Close et al., 2004; Janches et 
al., 2009). 

The Meteoroid Input Function (MIF) is a model developed to provide a global view of the 
diurnal and seasonal variations of the incoming meteoroid flux with the purpose of better 
understanding its relation to atmospheric phenomena in the MLT (Janches et al., 2006; 
Fentzke and Janches, 2008; Fentzke et al., 2009). The MIF assumes that the sporadic 
meteoroid population is distributed among the six known apparent sources: the North and 
South Apex, the Helion and Anti-Helion, and the North and South Toroidal (Jones and 
Brown, 1993; Taylor, 1997; Taylor and Efford, 1998). By assuming a global meteoroid mass 
input into the atmosphere (Ceplecha et al., 1998) and using Monte Carlo simulation 
techniques to assign an origin and speed to each modeled particle in agreement with the 
observed properties of the SMC apparent sources, the MIF predicts the overall meteoroid 
input to the MLT at any location on the Earth at any time of year and time of day. In order to 



validate the simulations with observations, the model also produces meteor detection rates of 
a particular HPLA system for a specified time and location. For this, considerations on 
sensitivity limitations of the specific HPLA instrument are utilized in order to characterize 
the portion of the incoming meteoroid population that is detectable by the specific radar 
(Close et al., 2005; Janches et al., 2008). An example of the large difference in distribution of 
meteor detections that can arise from using different observing instruments is shown by the 
meteor studies conducted with the 430 MHz Arecibo Observatory (AO) radar in Puerto Rico 
( ~ 1 8 0 N ) and the 440 MHz Poker Flat Incoherent Scatter Radar (PFISR) in Alaska (~67°N), 
in which both systems operate at a similar transmitting power and frequency yet have 
different aperture (Fentzke et al., 2009). The difference in radar aperture yields a considerable 
difference in detection sensitivity, and thus, there is a significant difference in the spectrum of 
meteoroid mass and velocity that each radar system is capable of observing. In addition, there 
is an extreme difference in geographic latitude between the two radar locations, which causes 
a difference in the portion of the SMC that each radar system observes. Nonetheless, the 
meteor detection rate of both HPLA systems and the seasonal and geographical variability of 
the incoming meteoroid flux were successfully captured by the MIF (Fentzke et al., 2009). 

In order to better assess the accuracy of the MIF and improve the model as needed, further 
comparison to meteor observations conducted by additional radar systems is necessary. In 
particular, utilizing systems with different transmitting power and frequency improves the 
confidence in the modeling of the physical considerations that determine radar sensitivity. 

For that purpose, this paper presents the simulation of meteor observations obtained with the 
46.5 MHz Middle and Upper atmosphere (MU) radar in Japan, offering the opportunity to 
validate the temporal variations modeled with the MIF at middle geographic latitude as well 
as the radar sensitivity considerations at the VHF frequency range. Lastly, the results 



presented here have the novel aspect that the MU radar system possesses the ability to obtain 
directional information from the detected meteors through the simultaneous use of several 
receiving channels, thus enabling the determination of the meteor radiant locations for the 
particle detections. This allows for the validation of the assumptions made regarding the 
sporadic meteor source and speed distributions used as inputs to the MIF. 

A description of the meteor experiment conducted at MU is presented in Section 2. The 
application of the MIF model and the characterization of the meteor detection sensitivity of 
the MU radar system are addressed in Section 3. Section 4 presents the results of the 
experiment by providing a comparison between the observed data and the model predictions. 
Lastly, the conclusions and final remarks are discussed in Section 5. 

2. MU RADAR METEOR OBSERVATIONS 

The MU radar is located in Shigaraki, Shiga Prefecture, Japan, at a geographic mid-latitude 
location (34.85°N, 136. 10°E) and has been utilized for over a decade for various studies 
involving meteor observations (Sato et ah, 2000; Fujiwara et ah, 2007; Kero et ah, 2011). The 
MU radar consists of 475 Cross-Yagi antennas arranged in a 103 m diameter circular array 
operating in the VHF band at 46.5 MHz and transmits a peak power of 1 MW. 

For the purpose of this work, meteor observations using the MU radar were conducted over 4 
distinct 24-hour periods near the equinoxes and solstices (December 18-19, 2008, and March 
20-21, June 18-19, September 24-25, 2009). These observation periods were specifically 
scheduled with the purpose of observing both the diurnal and seasonal variations in meteor 
activity at a mid-latitude geographic location. For the observations reported in this work, an 
uncoded radar pulse of 64 ps was transmitted with an interpulse period (IPP) of 1 .28 ms and 



a sampling frequency of 2 ps (i.e., 300 m range resolution). This scheme allowed for the 
probing of the altitude range between 79.5 and 139.2 km. In addition, the MU radar possesses 
25 receiving channels, enabling simultaneous reception with independent sub-sections of the 
antenna array and thus allowing for highly sensitive interferometric measurements 
(Hassenpflug et ah, 2008). For the case of meteor detections, this provides the capability to 
detennine the meteoroid’s vector position and velocity, which give orbital information for 
each detected event. For the observations presented here, however, only 4 receiving channels 
were utilized, but nevertheless interferometric capability was still maintained (Nishimura et 
ah, 2001). Figure 1 shows an example of a meteor event observed during the March 
campaign and summarizes the identification and analysis technique utilized. Figure la (upper 
left) shows the Range-Time-Intensity (RTI) plot from the signal recorded by one of the 
receiving channels (see Figure 1 in Nishimura et ah, 2001). The rectangular shape of the 
meteor head echo is due to the transmitted pulse shape (Sparks et ah, 2009). Figure lb (upper 
right) shows the Signal-to-Noise Ratio (SNR) of the same event. The overall ‘bell’ shape on 
the SNR curve is produced by the meteor traveling through the radar beam. It is important to 
note that the short-lived SNR enhancement occurring at approximately 0.03 s after the meteor 
is first detected is likely produced by the fast ablation of the alkali metals present in the 
meteoroid body (Vondrak et ah, 2008; Janches et ah, 2009). By cross-correlating the received 
raw voltages among the different channels, the phase differences between the received 
signals are obtained, as seen in Figure lc (middle left). This information is then utilized to 
calculate the elevation and azimuth of each IPP during which the meteor is detected 
(Nishimura et al., 2001; Chau and Woodman, 2004; Sparks et al., 2010), displayed in Figure 
Id (middle right). The line-of-sight velocity, which in this case is the vertical velocity, may 
be obtained from measuring either the rate of change of the meteor head echo altitude or the 
Doppler shift in the received signal. The absolute direction of the meteor is obtained from the 



measured range, azimuth, and elevation, as shown in Figure le and Figure If (lower left and 
right, respectively). This provides a complete set of parameters to allow orbit determination. 

The search and analysis of meteor events in the raw data follows the same methodology that 
has been applied to similar observations using the PFISR system (Sparks et ah, 2010), where 
the only difference is in the manner in which the Doppler radial velocity is estimated. For the 
440 MHz UHF operating frequency of PFISR, the Fast Fourier Transform (FFT) of the 
received signal in one IPP provides a sufficiently narrow spectrum to enable the accurate 
detennination of the returned signal Doppler shift (Sparks et ah, 2009). However, for the 46.5 
MHz operating frequency of MU, the spectrum is as broad as the radar bandwidth. The 
Doppler signature is then calculated by applying a pulse-to-pulse correlation function 
between consecutive IPPs (Hagen and Farley, 1973; Mathews, 1976; Janches et ah, 2000). 

3. MODEL DESCRIPTION AND APPLICATION 


3.1. The Meteoroid Input Function (MIF) 

The MIF is an empirical model developed for the purpose of providing a global description of 
the diurnal, seasonal, and geographical variability of the meteoroid input in the upper 
atmosphere. This knowledge is required to better assess the relationship between the 
incoming flux and the chemistry and composition of the MLT (Janches et ah, 2006; Gardner 
et al., 2011). The model, described in detail by Fentzke and Janches (2008), incorporates a 
global mass flux input into Earth’s atmosphere (Ceplecha et al., 1998) that is divided among 
the apparent sporadic meteoroid sources (Jones and Brown, 1993), each of which is 
represented as a Gaussian distribution of meteoroid speed and radiant. The North and South 
Apex sources are modeled collectively as approximately 33% of the meteoroids at 1 AU. The 



remaining two-thirds of incoming meteoroids is divided among the other apparent sources, 
with the Helion and Anti-Helion sources each contributing nearly 22%, and the North and 
South Toroidal each providing 1 1% of the total flux. The speed distribution for each source is 
based on radar observations from various studies, using both HPLA and meteor radar 
systems. The Apex source is modeled as bimodal with approximately 80% of the particle 
speeds centered at 55 km/s and 20% centered at about 17 km/s (Janches et al., 2003; Sulzer, 
2004; Janches and Chau, 2005). The remaining radiant sources (i.e., the Helion, Anti-Helion, 
and North and South Toroidal) are modeled as Gaussian with a peak at 30 km/s (Jones and 
Brown, 1993). Each of the apparent sources is modeled with the same particle mass 
distribution derived from the incoming mass-flux curve published by Ceplecha et al. (1998), 
which combines the results of several studies to estimate a global influx over the mass range 
of 10" 18 to 10 18 g and shows that the number of incoming meteoroids increases approximately 
exponentially with decreasing mass. The MIF models particles within the range of 10" 11 to 10" 
4 g; however, for the purpose of this work the results presented only include particles in the 
mass range of 0.1 to 100 pg. Particles in this range ablate and produce electrons independent 
of speed (Bronshten, 1983; Vondrak et al., 2008), but the ability to detect them is strongly 
dependent on the radar system utilized (see Section 3.2). Tables 1 and 2 summarize the 
modeled distributions, organized by sporadic source, of the meteoroids that are visible over 
the geographic latitude of the MU radar. The columns in Table 1 present the mean and 
standard deviation of the speed, ecliptic longitude (a), and ecliptic latitude ((3) distributions 
for the sources. The ecliptic longitude is the angle of rotation about the ecliptic normal 
measured from the Earth-Sun direction, and the ecliptic latitude is the angle of rotation out of 
the ecliptic plane (i.e., the Sun is located at k=0°, P=0°). Meanwhile, Table 2 provides the 
relative strength of contribution of each source, given as yearly and seasonal averages. 



Source 

Speed 

(km/s) 


H 

(deg) 

<*Spd 

(km/s) 

Ox 

(deg) 

<Tf 

(deg) 

Apex (Fast, Slow) 

55, 19 

270 

12 

6.6, 5.0 

19 

30 

Helion 

30 

349 

2 

6.6 

16 

16 

Anti-Helion 

30 

189 

3 

6.6 

18 

18 

North Toroidal 

30 

270 

58 

6.6 

19 

20 

South Toroidal 

30 

270 

— 

6.6 

16 

— 


Table 1: MIF modeled distributions of speed, ecliptic longitude, and ecliptic latitude for the 
meteoroids that are present above the geographic location of MU. 


Source 

Average 

Relative Contribution (%) 
Spring Summer Autumn 

Winter 

Apex 

29.8 

21 A 

29.3 

32.9 

29.0 

(Fast, Slow) 

(97, 3) 

(97, 3) 

(97, 3) 

(97, 3) 

(97, 3) 

Helion 

23.0 

25.0 

28.0 

21.6 

17.8 

Anti-Helion 

23.1 

26.3 

19.0 

19.5 

28.5 

North Toroidal 

20.0 

20.0 

19.6 

19.8 

20.5 

South Toroidal 

4.1 

1.2 

4.2 

6.2 

4.2 


Table 2: MIF modeled contributions of each sporadic source to the meteoroid population 

above the geographic location of MU. 


The MIF model provides the radiant location of each simulated particle, which is used to 
detennine the sporadic source from which the meteoroid originates. The meteoroid source is 
detennined by its proximity to the center locations of the six sporadic meteoroid apparent 
sources. The North and South Apex sources lie just above and below the direction of Earth’s 
velocity vector, with center points (X, P) at (270°, 20°) and (270°, -20°), respectively. 
Likewise, the North and South Toroidal sources are above and below the Apex sources and 
are centered at (270°, 60°) and (270°, -60°), respectively. The Helion direction is towards the 
Sun (0°, 0°), and the Anti-Helion is in the opposite direction (180°, 0°). In the MIF 
simulations, the apparent source center that is closest to the particle radiant is considered to 
be the source of the particle and is the basis for detennining the statistical distributions 
associated with each source presented in Tables 1 and 2. 


3.1.1. MIF Diurnal Rates of Meteoroid Areal Density 


The predicted rate of meteoroid counts per unit area, referred to here as the meteoroid areal 
density, at the MU radar latitude during the specified observation times is shown in Figure 2. 











It is important to note that the rates presented in the figure do not represent detected events 
but instead incoming particles in the mass range 0. 1-100 pg that form meteors at the specified 
geographic location (i.e., a meteor is formed but the MU radar may not detect it), assuming 
that the total meteoroid mass input is that given by the MIF and reported by Ceplecha et ah, 
(1998). From the panels in the figure, it is evident that there is a seasonal variability of the 
meteoroid areal density in both the total number and shape distribution of the diurnal rate. As 
expected, the peak of the daily meteoroid flux occurs during the autumnal equinox since this 
is when the Apex sources reach their highest elevation in the Northern Hemisphere sky. 
Conversely, the daily minimum occurs during the spring equinox when the highest elevation 
that the Apex sources reach during the day is lowest, with a large portion of the South Apex 
never reaching an elevation higher than the local horizon (Fentzke and Janches, 2008). In 
addition, when comparing the observation times of the winter and summer solstices, it can be 
seen that for the December period there is a sharper increase in the number of meteoroids in 
the early morning, while during the June period there is a steeper decrease in meteoroid 
activity in the late morning and early afternoon. This is due to the changing orientation of the 
Helion and Anti-Helion sporadic sources during the year with respect to the MU radar 
location, which is significant in determining the amount of activity around local midnight and 
noon. In December, the Anti-Helion radiant reaches higher elevation than the Helion in the 
Northern Hemisphere sky; whereas in June, the situation reverses, introducing seasonal 
variability around the midnight and noon periods (Janches et ah, 2006). In general however, 
at mid-northern latitude the simulated seasonal effects on the meteoroid input are not as 
pronounced as they are at higher latitudes (Janches et ah, 2006; Fentzke et al., 2009; Sparks 
et al., 2009). The peak rate for meteoroids in the 0. 1 to 100 pg mass range at the seasonal 
maximum (i.e., autumn) is estimated to be approximately 16 meteoroids/km 2 /min, compared 
to about 1 1 meteoroids/km 2 /min during the minimum activity period (i.e., spring). 



Finally, it is clear from Figure 2 that for each day there is a predicted peak in meteoroid 
activity at 6:00 am, local time, when the Apex reaches its highest point in the sky relative to 
the observer (i.e., the MU radar). The minimum in meteoroid areal density occurs at 6:00 pm 
for all observing periods due to the fact that most of the sporadic sources are below the local 
horizon at this time of day. Since the meteor observations conducted with MU are taken over 
the 4 different simulated 24-hour periods, it is expected that the features evident in the 
modeled diurnal variation of meteor activity will be prominent in the observed data. This 
analysis is presented in Section 4.2. 

3.1.2. MIF Modeled Distributions of Meteoroid Speed and Radiant 

In addition to the rate of meteoroid areal density, the simulated radiant and speed 
distributions for the incoming meteoroids are presented in Figures 3 and 4, respectively. The 
radiant maps given in Figure 3 show the modeled meteoroid radiant in tenns of ecliptic 
longitude (k) and latitude ((3) and represent the point in the sky that the meteoroids enter into 
a hyperbolic geocentric orbit (Jones and Brown, 1993). The radiant angles are defined as in 
Section 3.1, where the ecliptic longitude is the angle of rotation about the ecliptic normal 
measured from the Earth-Sun direction, and the ecliptic latitude is the angle of rotation out of 
the ecliptic plane (i.e., the Sun is located at k=0°, P=0°). The plots in Figure 3 are oriented 
such that the center point corresponds to the Apex direction (i.e., the direction of Earth’s 
velocity relative to the Sun, k=270°, (3=0°). The locations of the six sporadic meteoroid 
sources are also displayed in the figure as ellipses, with the coordinates as specified in 
Section 3.1. The North and South Apex sources lie just above and below the figure center 
point, respectively. Likewise, the North and South Toroidal sources are above and below the 
respective Apex sources. To the left of the Apex is the Helion direction, and the Anti-Helion 



is symmetrically opposite the Helion source about the Apex. Note that contributions from all 
of the sporadic sources are considered in the MIF and are clearly visible in Figure 3. 


It is important to note, once again, that the distributions presented in Figures 3 and 4, also 
summarized in Tables 1 and 2, are the results from the MIF previous to taking any 
consideration of the detection capability of the MU radar system. From these results it can be 
derived that, as expected, the Apex sources provide approximately 33% of the modeled 
meteoroids in the 0. 1 to 100 pg mass range, while the majority of the input is provided by the 
other sporadic enhancements. This is evident from the dominant peak in Figure 4 at 
approximately 30 km/s, as opposed to the secondary peak at 55 km/s, indicating that the 
majority of particles have the expected speed of the non- Apex meteoroids. 

It is also evident from the radiant maps in Figure 3 that, given a constant global meteoroid 
input, the relative contributions of the South Apex, Helion, Anti-Helion, and South Toroidal 
sources to the local MLT show a seasonal variability. For example, the Helion and Anti- 
Helion sources provide a similar contribution of meteoroids to the MLT at mid-northern 
latitude in March and September. However, during December and June, one source 
contributes more meteoroids than the other. This is consistent with the relative steepness in 
the increase and decrease of expected meteoroid activity in the observation dates for each 
season, as shown in Figure 2 and discussed in Section 3.1.1. Additionally, the radiant 
distributions in Figure 3 highlight the overall meteoroid activity for each observation period. 
In September, when meteor activity is at its peak in the Northern Hemisphere, there are clear 
contributions from each of the sporadic meteor sources, and there is a very strong 
contribution from the North Apex. On the remaining dates, it can be seen that some portions 
of the radiant map are cut off from appearing over MU, and the relative strength of the Apex 



sources is less than it is during September. The observation period in March demonstrates the 
inability to observe meteors with high southern ecliptic latitude radiant locations, and 
logically, this coincides with the time of year when meteoroid activity is at a minimum in the 
Northern Hemisphere. During this period, there is no contribution from the South Toroidal 
source and a reduced contribution from the South Apex source at mid-northern latitude. As in 
the case with the modeled variability of the diurnal rate of meteoroid areal density, these 
features are expected to be consistent with the MU observations, and the predicted and 
observed radiant maps will be compared to further validate the assumptions made by the MIF 
and the estimated input parameters used in the model. 

3.2. MU Radar Characterization 

When using instrumentation to conduct scientific experiments, it is necessary to characterize 
the biases and limitations introduced by the instruments being employed. For the case of 
meteor head echo observations with HPLA radars, it is inaccurate to generalize these system 
limitations over the entire radar class. Instead, each HPLA radar must be characterized and 
considered individually since the instruments are very different in comparison to each other 
(Janches et ah, 2008), and furthermore, the characteristics of each system are dependent upon 
several radar parameters (e.g., operating frequency, aperture, transmitted power, etc.) (Close 
et al, 2002). A similar approach must be used for any radar system, including meteor radars 
such as CMOR and AMOR (Brown et al., 2008). In this manuscript, we characterize the 
observational capability of the MU radar in terms of the minimum detectable radar cross- 
section (RCS). The meteor RCS can be detennined empirically from meteor measurements, 
assuming the meteor is a point target (Mathews et al., 1997; Close et al., 2002; Janches et al., 
2008; Fentzke et al., 2009), defined by the radar equation as 



RCS = 


64 Jt 3 P N p 4 
P t X 2 G 


SNR, 


(V 


In Equation 1 , p is the observed range to the meteor, G is the radar antenna gain, k is the 
signal wavelength, Pj is the radar transmitted power, and Pn is the noise power. The noise 
power is a function of the Boltzmann constant, k, the noise bandwidth, B, and the system 
temperature, T sys , and is given by PN=kT sys B- For the MU system, these values have been 
reported by Kero et ah, (2011) and are listed in Table 3. 


Parameter 

MU Radar System 

Px 

1 MW 

G 

2516 (peak) 

A 

6.45 m 

T sys 

12000 K 

B 

167 kHz 


Table 3: Operational parameters for the MU radar. 


In total, 12,430 meteors were identified in the raw data collected during the four observing 
periods. Figure 5 presents the speed versus altitude distribution for all of the collected data. 
The distribution shows a substantial peak in the observation count over an altitude range of 
95-105 km and at speeds between 50-60 km/s. This implies that, unlike the distributions 
shown in Figure 4, the vast majority of the detected meteors originate from the Apex apparent 
sources (Sekanina, 1976). Although this effect seems to be a typical result of HPLA 
observations (Close et al., 2000; Janches et ah, 2003; Westman et ah, 2004; Chau et al., 2007; 
Kero et al., 2011), the cause is different depending on the radar sensitivity (Janches et al., 


2008; Fentzke et al., 2009). 








The calculated RCS of the observed meteors by MU are shown in Figure 6. From the figure, 
a sensitivity threshold can be derived; that is, we can detennine a value for which meteors 
with smaller RCS will not be detected (Janches et al., 2008). From Figure 6, the threshold is 
estimated to be approximately -18 dBsm for the MU radar. This detection limit is 
significantly higher (i.e. the MU radar is less sensitive) than the thresholds that have been 
detennined on other HPLA radars that are operating at higher frequencies (Close et ah, 2007; 
Janches et ah, 2008; Fentzke et ah, 2009). Additionally, Kero et ah (2011) reports much 
lower RCS values for MU observations conducted during a separate experiment, which are 
due to the different operating conditions of the MU radar system prevalent in each study. This 
highlights the necessity to assess the sensitivity characteristics of each radar system 
individually and for the specific operating parameters utilized; generalizing the radar 
characteristics to other systems and observational studies can lead to inaccurate results. 

Once the radar RCS threshold is determined, the properties of the meteoroids that are 
observable by the MU system are estimated by applying a scattering model that considers the 
meteor head echo as a spherical target with Gaussian plasma density distribution (Close et al., 
2004; 2005; 2007). Assuming that the incident signal wavelength is larger than the head 
plasma radius (i.e., the radar frequency is at UHF or below) and the incoming meteoroid is a 
solid sphere, the meteor RCS is calculated as a function of particle mass, speed, and ablation 
altitude for a given radar operating frequency, described in detail by Close et al. (2005; 

2007). Shown in Figure 7 are the modeled results when applied to the MU radar system. In 
order to highlight the differences among HPLA radars, Figure 7 also shows the results for 
three additional HPLA systems used for meteor observations: the ARPA Long-range 
Tracking And Instrumentation Radar (ALT AIR; 160 MHz), Arecibo (430 MHz), and PFISR 
(440 MHz). The vertical line depicts the sensitivity threshold of the radar, which is 



determined for MU using the observed data presented in Figure 6 and for the other radar 
systems using similar methodology presented in previous studies (Close et ah, 2007; Janches 
et ah, 2007; Fentzke et al., 2009). According to the model, particles with RCS values that are 
to the right of the sensitivity threshold in Figure 7 will be detected by the respective radar, 
while meteoroids with RCS to the left will not produce a strong enough signal to be detected 
(Janches et ah, 2008). In comparing the sensitivity limits of the different instruments, some 
important conclusions can be made. First, the RCS of a meteoroid with given mass and speed 
is highly dependent on the radar operating frequency. This dependence scales approximately 
as f 2 for meteor head echoes and f 6 for meteor non-specular trails, where f is the radar 
operating frequency (Close et al., 2008). Second, the sensitivity threshold varies greatly 
among the different HPLA systems. The MU radar, operating at 46.5 MHz, has a threshold at 
approximately -18 dBsm, which is much higher than the threshold of the 430 MHz Arecibo 
system, for example, at approximately -90 dBsm (Janches et al., 2007; Fentzke et al., 2009). 
Thus, it is clear from Figure 7 that the HPLA radar systems observe significantly different 
portions of the sporadic meteoroid population, and the unique biases of the particular 
instrument being utilized must be characterized in order to appropriately assess the properties 
of the incoming particles. Lastly, it is to be noted that this RCS modeling technique has been 
utilized in previous studies to produce consistent estimates of meteoroid mass from 
simultaneous meteor head echo observations on VHF and UHF with the ALT AIR system 
(Close et al., 2005; 2007). Additionally, the application of the RCS sensitivity threshold to 
the MIF model has been shown to accurately predict meteor observations conducted at 
Arecibo and PFISR (Fentzke et al., 2009). 


The RCS modeled results suggest that the MU radar efficiently observes meteoroids with 
masses on the order of 1 pg if they are moving with speeds greater than about 60 km/s. For 



larger particles (~10 pg), detectable RCS values will be produced at lower speeds (~30 km/s). 
For masses greater than 10 pg, all particle speeds should produce meteors with RCS 
observable by MU. However, the occurrence rate for such objects will be much lower than 
that for the less massive particles (Ceplecha et ah, 1998). This is supported by the observed 
meteor speed distribution shown in Figure 5, which indicates that the majority of meteors 
detected by MU have speeds within the range 50-60 km/s. This situation is much different 
when considering other radars, however. The Arecibo 430 MHz radar system, for example, is 
able to effectively observe meteoroids as small as 0.1 pg over almost the entire incoming 
speed range and can observe even smaller masses if the meteor speed is higher. 


Mass 


Minimum Speed (km/s) 


Gog log) 

MU 

ALTAIR 

Arecibo 

PFISR 

-7 

80 

40 

25 

— 

-6 

60 

25 

15 

25 

-5 

25 

15 

5 

15 

-4 

10 

0 

0 

0 

-3 

10 

0 

0 

0 


Table 4: Minimum meteoroid speed required for radar detection as a function of meteoroid 

mass for several HPLA radar systems 

The speed thresholds of observable meteors for the MU, ALTAIR, Arecibo, and PFISR 
systems are listed in Table 4 for the range of meteoroid masses that is included in the MIF 
model for this study. Note that many of the meteoroid sizes included in the MIF simulation 
presented in Section 3.1 are only observable by the MU radar if the particles are traveling 
with speed greater than approximately 50 km/s. Moreover, particles smaller than 0. 1 pg 
regardless of speed are completely undetectable by the MU radar when using this particular 
wavefonn. Table 4 highlights that, once again, one must be cautious about the misleading 
interpretations that may occur when observation results from a particular HPLA system are 
taken out of context or generalized to other HPLA radar systems (Hunt et ah, 2004; von 











Zahn, 2005). While the analysis of the MU radar shows that a 1 jag meteoroid will not be 
detected if its incoming speed is slower than 60 km/s, the same methodology shows that a 
particle of this mass will be observed at much lower speeds by separate, more sensitive 
HPLA instruments (Figure 7; Janches et ah, 2008; Fentzke et ah, 2009). 

Understanding the instrument dependent selection bias is crucial to estimating the total input 
of meteoric material in the atmosphere from these observations. For example, Nesvorny et al. 
(2010; 2011) report results from a recently developed dynamical model of the Zodiacal Dust 
Cloud (ZDC). The model considers the orbital properties of comets and asteroids from the 
dynamical evolution of dust particles after ejection, and its results are constrained by Infrared 
Astronomical Satellite observations of the ZDC. The authors suggest that 85-95% of the 
continuous input of sporadic meteoroids comes from Jupiter Family Comets, with the 
remainder originating from the asteroid belt and Oort cloud comets. Additionally, this portion 
of the continuous input has mass in the range 1 to 10 pg and a mean speed of about 14 km/s, 
providing a global mass input of roughly 42x1 0 3 kg/d. Furthermore, the authors suggest that, 
due to their slow speed, most of the incoming particles will not produce significant ioni z ation 
to be detected by radar systems, and thus, this large flux occurs largely undetected by ground- 
based instruments. A qualitative attempt to validate this conclusion using meteor radar 
observations by the Canadian Meteor Orbit Radar (CMOR; Webster et al., 2004) and the 
Advanced Meteor Orbit Radar (AMOR; Baggaley et al., 1994) is reported in Nesvorny et al. 
(2011). Additionally, the results presented in Figure 7 and Table 4 demonstrate that the MU 
system would not be able to detect the majority of particles entering the atmosphere at slow 
speeds with mass between 1 to 10 pg, and if indeed the population proposed by Nesvorny et 
al. (2010; 2011) exists, it would remain mostly undetected by this particular radar system. 
However, when considering the results given in Figure 7 and Table 4, it is evident that such a 



meteoroid population would be detected by the Arecibo radar, and the ALTAIR and PFISR 
systems would be able to detect a significant portion of it. This implies that the speed 
distribution measured by the most sensitive HPLA radar systems would have to be dominated 
by slow particles, and yet this has not been observed to be the case. 

4. COMPARISON OF THE OBSERVATIONS TO MODEL SIMULATION 

Utilizing the MIF results to reproduce the meteor observations of HPLA radar systems has 
shown great agreement in previous work. This requires the implementation of methods to 
characterize the detection characteristics of the radar systems employed. Specifically, meteor 
observations made with Arecibo in Puerto Rico and PFISR in Alaska have been accurately 
modeled using the MIF together with the empirically detennined RCS threshold of each radar 
system (Fentzke et ah, 2009). Meteoroids with masses and speeds that coincide with a meteor 
RCS less than the sensitivity threshold of each instrument (Section 3.2) are considered 
undetectable and removed from the MIF output since the model predicts the total meteoroid 
presence without any consideration of specific observational biases (Section 3.1). The 
resulting distribution from applying this method represents the simulated meteors detected in 
each experiment. With both Arecibo and PFISR, the predicted meteor detection rate and 
radial speed distribution strongly agreed with the observed results. Furthermore, the 
simulations also showed strong agreements in the diurnal, seasonal, and geographical 
variability of these quantities. Additionally, the results show that even though the particular 
method of determining the RCS sensitivity limit of two unique HPLA radar systems 
represents a parameterization of a complex process that does not yet consider some effects 
(e.g. radar beam pattern or the latest ablation and electron production models), it nevertheless 
produces accurate results in predicting the overall observed meteor rates and speed 
distributions. In this section, the same methodology is applied with the MIF simulations 



presented in Section 3. 1 in conjunction with the RCS sensitivity results of the MU radar 
described in Section 3.2 in order to compare the model-predicted meteoroid properties with 
the meteor observations conducted at MU. 

4.1. Meteoroid Entry Angle 

A technical advantage for meteor observation that is available with the MU radar is its 
capability to perform interferometry, enabling the detennination of the 3 -dimensional 
position and velocity, and thus the orbit, of a meteoroid. Likewise, the zenith angle of the 
meteoroid trajectory can be obtained from the components of the velocity vector and is 
referred to here as the meteoroid entry angle. The observed entry angle distributions for the 
MU observation periods are shown in Figure 8. In the figure, an entry angle of 0° 
corresponds to a trajectory that is aligned with the local vertical (i.e. the meteoroid is 
travelling straight downward), while 90° corresponds to a horizontal velocity vector. A clear 
seasonal effect is evident in these distributions. As expected from previous work, the entry 
angle distributions are similar during the solstices, while during the equinox dates the angles 
are skewed in opposite directions from the solstice distribution. At the autumnal equinox, the 
majority of meteoroid trajectories are closer to vertical, while during the spring equinox the 
trajectories are more horizontal. These results are in agreement with previous observations, 
and the seasonal effects have been observed to be more pronounced at higher geographic 
latitudes (Janches et ah, 2006; Fentzke et ah, 2009; Sparks et ah, 2010). This effect is due to 
the relative position of the apparent sources, in particular the Apex, with respect to the 
observer (i.e. the MU radar). Characterizing this distribution of meteoroid entry angle enables 
the construction of a more realistic picture of how the incoming meteoroid mass is deposited 
in the upper atmosphere and how it impacts the composition and chemistry of the mesosphere 
(Gardner et ah, 2011). Additionally, the distributions presented in Figure 8 allow for the 



estimation of a range of entry angles within which the incoming meteoroids detected by MU 
appear to occur. A similar methodology has been used in previous work employing MIF 
simulations, where meteoroids with radiant locations below 20° elevation from the local 
horizon appeared not to be observed with the respective radar systems (Janches et ah, 2006; 
Fentzke and Janches, 2008; Fentzke et ah, 2009). However, for the MU radar, the entry angle 
of each detected meteoroid can be calculated, and the range of observed entry angles can be 
empirically determined with more confidence. Additionally, the entry angle range can be 
estimated individually for each observation period to provide a more accurate description of 
the observation capabilities of the MU radar. In order to determine the entry angle range, the 
points where the number of observations with a particular entry angle is half of the number of 
observations of the most observed entry angle (i.e., the entry angle values in Figure 8 where 
the distribution crosses 0.5 on the vertical axis) are taken as the limits to the range. These 
angle ranges for the MU observation periods were detennined to be: 33°-65° for March; 25°- 
60° for June; 17°-60° for September; and 20°-60° for December. The characterization of the 
MU radar entry angle range is necessary because meteoroids outside of this range appear in 
the MIF simulations, and these particles must be removed from the model in order to 
accurately estimate properties of the meteors detected with the MU radar system, such as the 
areal density rate, speed distribution, and radiant positions. 

4.2. Meteor Areal Density Rate 

In Section 3.2 the MU meteor observations were utilized to characterize the sensitivity of the 
radar system. For the purpose of this work, the mass-speed threshold presented in Figure 7 
and summarized in Table 4 is used to determine which of the modeled meteoroids will be 
detected by the MU radar. As described in Section 3.2, meteoroids with mass and speed 
characteristics that are to the right of the RCS thresholds in Figure 7 are considered 



observable by that particular radar system, which enables the detennination of the minimum 
meteoroid speed that a particle with a given mass must have in order to be detected. All 
modeled meteoroids with a mass and speed that is below the detection threshold are removed 
from the MIF simulation of Section 3.1. In addition, meteoroids originating from a radiant 
that has an elevation angle outside the range of entry angles that are observed to produce the 
majority of the detections, as shown in Figure 8, are neglected from the simulated input. We 
note once again, that one must not generalize a particular mass-speed threshold to all HPLA 
radar systems. The sensitivity characteristics are only applicable to the specific radar 
instrument for which they were determined. 

Comparisons between observed and predicted meteor areal density rates (i.e., the number of 
detected meteoroids per cross-sectional area per unit time) for the four observation periods 
are presented in Figure 9. It can be seen from the figure that there is very good agreement 
between the observed and predicted meteor rates. Note that these numbers are not normalized 
and represent actual detection rates. Moreover, these results indicate that even though the 
simulated detection threshold methodology is somewhat simplistic, it accurately reproduces 
the observations, which peak between 0.2-0. 3 meteors/km 2 /min for each of the observation 
periods. Further development of the sensitivity threshold analysis is ongoing and will include 
head echo electron density profiles derived from differential ablation models (Vondrak et al., 
2008; Janches et al., 2009) as well as a more accurate representation of the radar beam shape 
(Dyrud and Janches, 2008). 


As expected, for each observing period, the observed and predicted rates of areal density peak 
around 6:00 am, when the Apex reaches its highest elevation. The seasonal minimum rate of 
meteor areal density occurs in March (i.e. spring), resulting in a significantly lower meteor 



rate than the other observing periods. This is also predicted by the model and is consistent 
with the overall trend seen in the modeled meteoroid input at the MU system’s geographical 
location, illustrated in Figure 2. Note that the maximum rate of observed meteors occurs in 
December and not in September, which is unlike the MIF input that predicts the maximum 
rate of activity during the autumn. A similar trend was reported in Arecibo observations 
(Janches et ah, 2006). One characteristic of the MIF is that it assumes a uniform meteor input 
throughout the year, and these results indicate that this may not be the case. For some times 
of year, the overall sporadic meteoroid background input may be larger than at other times. 
Regardless of these unexpected trends, in general the MIF results shown in Figure 9 
accurately predict the seasonal variations observed by the MU radar when also accounting for 
the sensitivity characteristics of the system. This highlights the ability of the model to 
detennine the input of meteors and the corresponding radar-dependent portion of the 
meteoroid population that is detectable. These detection biases and limitations of the 
particular radar system can introduce additional diurnal and seasonal variations to the 
observed data, which must be characterized accurately in order to properly model the impact 
of the sporadic meteoroid population in the upper atmosphere. For example in Figure 9, the 
December and June observations show skewed peaks in the meteor activity around the time 
when the elevation of the Apex is highest (i.e. 6:00 am local time). There are more meteors in 
the early morning in December and in the late morning in June. This is also captured in the 
MIF simulations (Figure 2) and is due to the orientation of the radar location with respect to 
the Helion and Anti-Helion sporadic sources at these times. This effect becomes more 
pronounced when the detection characteristics of the MU radar are applied to the MIF results. 
This causes the skewed peaks in the predicted detection rates evident in Figure 9 for 
December and June as well as the double hump structure in March and September. These 
same features have been observed at lower latitude by Arecibo (Janches et ah, 2006). In 



December, the maximum elevation of the Anti-Helion source is higher than that of the 
Helion, and thus, more Helion originating meteoroids will be below the elevation threshold 
(Section 4.1) and rejected by the model. This causes a dip in the late morning detections in 
December. The reverse is true in June, causing the dip to occur in the early morning. The 
double hump in the detections during March and September are due to the position of the 
Apex sources in the sky. In September, the Apex is better visible in the Northern Hemisphere, 
and so the effect is more pronounced, as more of the Apex meteoroids are above the 
maximum elevation that tends to be observed, as represented in Figure 8. 

The major difference between the model predictions and the radar observations is the 
overestimate by the model of the meteor areal density rate in the late morning and early 
afternoon for three out of the four observing periods. During these three periods, the modeled 
downward slope in areal density rate shown in Figure 9 appears to be consistent with the 
downward slope in the observed rate, but the observations start to decrease approximately an 
hour before the model predictions. Further analysis is required to determine the cause of this 
difference between the MIF model results and the observations. Nonetheless, these results 
show that the application of the MIF model together with the parameterization of the 
particular radar system used for the observations accurately captures the seasonal and diurnal 
variations in detection rate at mid-latitudes as it has for equatorial and polar latitudes using 
Arecibo and PFISR, respectively (Fentzke et ah, 2009). In addition, these results demonstrate 
that the model can be applied to much lower radar operating frequencies than have been used 
in previous work. 



4.3. Meteor Speed Distribution 


An important property needed to understand the impact of sporadic meteoroids in the MLT is 
the speed distribution of the incoming particles. It is crucial to understand the differences 
between the true input speed distribution and that which is observed by a particular radar 
instrument. For the case of head echo observations using HPLA radar systems, this has been 
a fiercely debated topic, mainly due to the tendency of generalizing the results of one 
particular system to the entire radar class (Close et ah, 2007; Janches et ah, 2008). The input 
speed distribution considered in the MIF model is discussed in Section 3.1.2 and shown in 
Figure 4. The observed results by the MU radar as well as the instrument-specific predicted 
distributions, which account for the radar detection biases, are shown in Figure 10. Since the 
MU radar observations utilize interferometry, the speeds given in Figure 10 are the 
magnitudes of the absolute velocity vectors of the detected meteoroids and not the line-of- 
sight components as presented in previous studies (Fentzke and Janches, 2008; Fentzke et ah, 
2009). From Figure 10, it is evident that there is remarkable agreement obtained between the 
observed and predicted distributions. During each of the observing campaigns, the 
observations are characterized by a dominant peak centered at approximately 55 km/s, which 
is also well reproduced by the model. This indicates that the observations are dominated by 
particles originating from the Apex sources and is consistent with the observed maximum in 
diurnal rate of meteor areal density during the time when the Apex is at its highest elevation, 
as shown in Figure 9. 


MU SPECIFIC MIF MODEL PREDICTIONS 


Source 

Speed 

X 

ii 

&Spd 

Ox 

Of 

(km/s) 

(deg) 

(deg) 

(km/s) 

(deg) 

(deg) 

Apex (Fast, Slow) 

57, 25 

270 

17 

6.2, 6.3 

19 

25 

Helion 

34 

350 

5 

7.0 

16 

14 

Anti-Helion 

34 

189 

5 

6.9 

18 

16 

North Toroidal 

35 

270 

58 

6.9 

20 

20 

South Toroidal 

35 

266 

— 

7.0 

14 

— 


MU METEOR OBSERVATIONS 



Source 

Speed 

X 

n 

&Spd 

Ox 

Of 

(km/s) 

(deg) 

(deg) 

(km/s) 

(deg) 

(deg) 

Apex (Fast, Slow) 

59, 29 

270 

6 

10.2, 5.3 

16 

17 



Helion 

37 

346 

4 

19.5 

8 

14 

Anti-Helion 

35 

194 

3 

15.3 

1 

13 

North Toroidal 

48 

270 

54 

16.9 

18 

10 

South Toroidal 

54 

266 

-48 

24.4 

16 

6 


Table 5: Modeled and observed distributions of speed, ecliptic longitude, and ecliptic latitude 
for the meteoroids detected by the MU radar system. 


Table 5 displays the mean and standard deviation in observed speed of the incoming 
meteoroids originating from each source. The determination of the meteoroid source for both 
the observed and modeled meteors was conducted in the same manner as the source 
detennination of the MIF simulated particles described in Section 3.1. The table shows that 
there is good agreement in the predicted and observed speed distributions of meteoroids from 
the Apex sources, where the predicted mean speeds are within 3.5% and 15% of the observed 
speeds for the fast and slow components, respectively. There is a significant difference in the 
standard deviation of the fast component, however, where the predicted value is almost 50% 
lower than the observed value. This is also the case with the observed standard deviation in 
speed of the other sporadic sources. The predicted deviations are 50-70% less than the 
standard deviations observed by the MU radar system. Moreover, some unexpected results 
are shown in the table between the predicted and observed mean speeds for these sources. 
While the contributions of the Helion and Anti-Helion sources are shown to be significantly 
overestimated by the MIF simulations, the estimated mean speed of the meteoroids 
originating from these two sources agrees well with the observed values. However, this same 
agreement is not evident between the predicted and observed values for the Toroidal sources. 
As in the case of the Helion and Anti-Helion sources, the MIF simulations predict that the 
mean speed of incoming meteoroids from the Toroidal sources is approximately 35 km/s; 
whereas the observations show a mean speed of 48 km/s and 54 km/s for the North Toroidal 
and South Toroidal, respectively. However, since the Toroidal sources are relatively close to 
the Apex sources, it is likely that this result is due to a contamination on the Toroidal sample 





by particles originally belonging to the Apex sources (i.e., these particles originate from a 
radiant location where the distributions of the Apex and Toroidal sources overlap). These 
results suggest that better orbital constraints, increased statistics, and further analysis and 
modeling of the speed characteristics of the sporadic sources are required, particularly over 
the meteoroid mass range of 0. 1 to 100 pg that the MU system observes. 

As shown in Table 4, the minimum speed for a 1 pg meteoroid to be detected by the MU 
radar is 60 km/s, but for slightly larger masses the minimum speed threshold decreases 
rapidly. Since the mean speed of the Apex source occurs near this threshold speed and the 
modeled number of meteoroids increases exponentially with decreasing mass (Ceplecha et 
ah, 1998), this means that the observations (and predictions) will be dominated by the 
smaller, faster particles. Some of the differences between the shapes of the observed and 
predicted distributions are likely introduced by the fact that small errors in either the mass 
distribution or the radar sensitivity model can lead to fairly large changes in predicted 
meteoroid areal density and speed. Finally, there also appears to be a small seasonal variation 
in the peak location and width of the speed distribution that will need to be incorporated into 
the model to improve the predictions. The seasonal variation is even more prevalent at the 
lower speeds due to the fact that the relative strength of the different sources is dependent on 
time of year. Uncertainties in the mass distribution and radar threshold at these lower speeds 
likely explain the differences between the modeled and observed distribution. Furthennore, 
there may also be errors in the modeled speed distribution of the non- Apex sources that 
account for these observed differences, particularly with the Helion source in which the 
predicted areal density rate at times associated with this source also showed slight 
differences. Nevertheless, these errors appear to be mainly higher order, and the general 
shape of the predicted speed distribution agrees with the observations very well and is 



sufficiently accurate for practical purposes. One final note is that such differences were not 
evident in previous observations conducted with the Arecibo radar, which detects much 
smaller particles than MU (Janches et ah, 2006; Fentzke and Janches, 2008). In particular, the 
AO system is able to detect particles with mass as low as 0. 1 jig traveling at slow speeds. 

This indicates that the sporadic source characteristics are likely dependent on meteoroid 
mass. 

4.4. Meteor Radiant Distribution 

The advantage of perfonning interferometry measurements with the MU radar is that it 
allows for the determination of the meteoroid orbital properties and radiant location (i.e., the 
point in the sky in which the meteoroid appears to originate; Kero et ah, 2011). Table 6 
summarizes the predicted and observed contributions of the sporadic meteoroid sources to the 
MU meteor detections and includes the radar-specific sensitivity characteristics, as discussed 
in Section 3.2. Once again, the determination of the observed and modeled meteoroid 
sources used the same methodology as described in Section 3.1. The table shows that the MIF 
model obtains strong agreement with the observations when the radar detection 
characteristics are taken into account. The strong effects introduced by the sensitivity of the 
MU system are further evident when comparing these results to the modeled input source 
contributions listed in Table 2. As discussed in Section 3. 1, for meteoroids with masses in the 
0.1 to 100 pg range, the Apex sources collectively contribute approximately 33% of the MIF 
model input, while the remaining 67% originates from the other sporadic sources, mostly 
from the Helion and Anti-Helion. However, as summarized in Table 6, the Apex contributes 
approximately 77% and 87% of the predicted and observed meteors, respectively, that are 
detected by the MU radar. While there is a substantial difference between these two values, 
both the predicted and observed contributions indicate that the vast majority of the MU- 



detected meteors originate from the Apex sources, with the remaining sporadic sources 
representing a minority of the detections. The results in Table 6 show that the modeled 
contribution of the Helion and Anti-Helion sources is overestimated significantly, with 
approximately 20% of the detections predicted to come from these two sources collectively 
and only about 5% of the observed meteors attributed to them. Furthennore, while the 
modeled and observed contributions of the North Toroidal source agree strongly, this source 
is the secondary contributor to the meteors observed by MU yet is predicted to be only the 
fourth strongest sporadic source by the MIF simulations. 


MU SPECIFIC MIF MODEL PREDICTIONS 


Source 

Average 

Relative Contribution (%) 
Spring Summer Autumn 

Winter 

Apex 

76.6 

68.3 

77.3 

79.0 

78.7 

(Fast, Slow) 

(99,1) 

(99,1) 

(99,1) 

(99,1) 

(99,1) 

Helion 

8.5 

13.5 

12.2 

7.6 

4.1 

Anti-Helion 

9.0 

12.7 

3.7 

7.7 

11.3 

North Toroidal 

5.6 

5.2 

6.5 

5.3 

5.5 

South Toroidal 

0.3 

0.3 

0.3 

0.4 

0.3 

Source 

MU METEOR OBSERVATIONS 

Relative Contribution (%) 
Average Spring Summer Autumn 

Winter 

Apex 

86.8 

80.8 

87.8 

86.2 

90.0 

(Fast, Slow) 

(98,2) 

(98,2) 

(98,2) 

(98,2) 

(98,2) 

Helion 

2.8 

4.9 

3.6 

1.8 

2.1 

Anti-Helion 

3.0 

5.9 

2.7 

2.6 

2.0 

North Toroidal 

6.6 

8.3 

5.4 

7.6 

5.3 

South Toroidal 

0.8 

0.0 

0.5 

1.8 

0.5 


Table 6: Modeled and observed contributions of the sporadic sources to the meteoroid 
population detected by the MU radar system. 


Figures 1 1 and 12 display the predicted and observed radiant distributions, respectively, for 
the meteors detected by the MU radar during the four observation periods, which are also 
summarized in Table 6. The coordinates in these figures as well as the location of the 
apparent sources are the same as in Figure 3 and described in Section 3.1.2. While it is more 
difficult to visually compare the meteor radiant distributions than it is the distributions of 
meteoroid areal density rate and speed due to the fact that the predicted and observed radiant 
distributions are plotted separately and span an additional dimension, some important 









statements can nevertheless be made. In both the modeled and observed distributions, it can 
be seen that, as expected, the majority of the meteoroids originate from the Apex sources, 
particularly the North Apex. This is consistent with past studies as well as the observed 
variations in meteor areal density and speed (Sections 4.2 and 4.3) and agrees with the 
relative strength of contribution of the Apex sources given in Table 6. However, the predicted 
distributions of the Helion and Anti-Helion sources show stronger contributions than what 
was observed. This may suggest that the modeled strength and deviation of these sources is 
too high, or as discussed in Section 4.3, the mass and speed distributions of these sources 
may need to be reevaluated. Finally, the seasonal variation in the predicted meteor radiant 
distribution is clearly evident in Figure 11, where the portions of the ecliptic coordinates that 
are unobservable from the MU radar location are easily identified. While not as clearly 
visible in the MU meteor observations, the same seasonal variation is seen in the observed 
radiant distributions in Figure 12. The portion of the radiant map that is cut off in the 
predicted distribution is also free of observed radiant locations for each observation date. This 
shows the ability of the MIF to predict the general origin of meteoroids that are present above 
specific locations on the Earth at different times of year. 

Similar to the discussion in Section 4.3, the results presented in Table 5 allow for some 
important comparisons to be made between the predicted and observed spatial distributions of 
each of the sporadic sources. The table provides the mean and standard deviation of the 
ecliptic latitude and longitude (defined in Section 3.1) of each source predicted by the MIF 
simulations and observed during the MU meteor observation periods. For each source, the 
observed mean ecliptic longitude (X) agrees well with the predicted values. The largest 
difference is with the Helion and Anti-Helion sources, where the observed mean ecliptic 
longitude is approximately 5° closer to the Apex than that which is given by the MIF results. 



Likewise, the observed standard deviation in ecliptic longitude shows strong agreement with 
the predicted values for each of the sporadic sources except for the Helion and Anti-Helion. 
For these two sources, the observed distribution shows very little variation in ecliptic 
longitude, while the model suggests a standard deviation of 15-18°. For ecliptic latitude, [3, 
the results are much different. The predictions and observations show excellent agreement in 
both the mean and standard deviation of ecliptic latitude for the Helion and Anti-Helion 
sources. However, for the remaining sporadic sources, the MIF simulations tend to place the 
expected ecliptic latitude further north than what is observed by the MU radar. The model 
also predicts a much larger standard deviation in ecliptic latitude than what was observed. 
Clearly, more observations are required to improve the statistics in the weaker sporadic 
sources, and additionally, further analysis is needed to shed light into the differences between 
observed and predicted meteoroid speeds and radiant locations. 

5. CONCLUSIONS 

In this paper, the analysis of meteor head echo observations conducted with the 46.5 MHz 
MU radar in Japan is presented. The observation periods were specifically scheduled with the 
purpose of observing both the diurnal and seasonal variations in meteor activity at a mid- 
latitude geographic location over four distinct 24-hour periods near the equinoxes and 
solstices (December 18-19, 2008; and March 20-21, June 18-19, September 24-25, 2009). 

The observed meteor properties were derived with the MIF model developed by Fentzke and 
Janches (2008) in a similar manner utilized in previous meteor studies performed with the 
Arecibo and PFISR systems (Fentzke et al., 2009). However, unlike previous studies, the 
multiple receivers of the MU radar enabled interferometry measurements, allowing the 
derivation, comparison, and validation of additional meteor flux characteristics. In order to 
accurately compare the observations with the MIF, the sensitivity of the MU radar was 



parameterized in order to determine the portion of the incoming meteoroid population that is 
detectable by the radar system in terms of particle mass and speed. A comparison of the MU 
radar parameterization with that of the HPLA radar systems ALT AIR, Arecibo, and PFISR 
was also presented. This direct comparison of the sensitivity characteristics of the different 
instruments demonstrates that each radar system detects a very distinct mass-speed portion of 
the incoming meteoroid flux. This is a very important result in light of a recently developed 
model of the Zodiacal Dust Cloud that predicts a large amount of meteoric mass in the 
atmosphere that remains undetected by ground-based systems due to the slow speeds of the 
incoming particles (Nesvorny et ah, 2010; 2011). The results presented here show that 
although this model would be consistent with the sensitivity characteristics of the MU radar, 
more sensitive HPLA systems would be able to detect some portion of these slow particles. 

Once the parameterization of the MU system’s sensitivity is applied to the MIF, the modeled 
results show strong agreement with the observed seasonal and diurnal variations of meteor 
areal density over the course of the observation periods. This result validates the accuracy of 
the model at mid-northern latitudes. In addition, given the MU system’s VHF operating 
frequency, the methodology used in this paper can, as a result, be applied to lower 
frequencies than previously demonstrated. 

Besides the diurnal rates, the MIF accurately reproduces the meteor speed distribution 
observed by the MU radar. This supports the modeled input characteristics of the MIF for the 
range of masses that are observable by the radar. As expected, the results show that the 
majority of the meteor detections originate from the Apex sources and have masses as low as 
1 pg for high speed particles (~60 km/s). While there is strong agreement in the distribution 
at faster speeds, there is less agreement at the slower meteor speeds. This is probably due to 



the fact that at the slow speeds the overall meteor counts are lower, and thus, the errors are 
more pronounced. However, it may also suggest that further analysis is required in the 
modeling of the meteoroid mass and speed distributions associated with the non-Apex 
sources and their relative strengths. 

The MIF has also been shown to accurately predict the observed meteor radiant distributions. 
The seasonal variation in the origin of the observed particles is correctly modeled, and the 
relative strengths of the sporadic sources appear accurate. However, the distribution within 
each source seems to require more analysis as the MIF predicts a wider standard deviation in 
the radiant locations than what was observed. Nevertheless, the comparison of the model 
predictions to the MU observations greatly supports the validity of the MIF and its ability to 
be applied to different locations on the Earth at any time. 

While the MIF has been shown to be very effective in predicting the properties of the 
meteoroid population that are observable with HPLA radars, it is important to note additional 
studies that are needed to enhance the strength of the model. First, the HPLA radars that have 
been used in MIF studies have all employed zenith pointing antennas. To improve the 
confidence in the MIF, it would be useful to compare its perfonnance in predicting the 
meteor detections of radar systems that are observing different regions of the sky, such as 
those by the ALT AIR system, which has been used for meteor observations and has a 
steerable radar antenna (Close et ah, 2000; Hunt et ah, 2004). This would allow for better 
refinement of the sporadic source distributions. Second, the radar observations that have been 
compared to the MIF have been primarily from northern geographic locations. Although 
Janches et al. (2006) show agreement between an earlier version of the MIF and observations 
with the Jicamarca Radio Observatory (JRO) in Peru, it would be valuable to use the current 



version of the model and compare it with additional meteor observations from the Southern 
Hemisphere. These studies can be performed not only with JRO but also with the Southern 
Argentine Agile MEteor Radar (SAAMER) in Argentina (Fritts et ah, 2010), which is located 
at high-southern latitude. Such instruments would provide insight into differences in meteor 
activity over the Northern and Southern Hemispheres, particularly in regard to any 
asymmetry in strength between the North and South Toroidal sporadic sources. In addition, 
potential studies with SAAMER provide the opportunity to apply the MIF model to meteor 
radar observations of specular trails, which will greatly increase the confidence in the 
accuracy of the model, particularly for larger and slower meteoroid masses originating from 
the non-Apex sources. 
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8. FIGURE CAPTIONS 


Figure 1 : Example meteor event detected by the MU radar system during the March 
observations - a) Range-Time-Intensity (RTI) plot taken from one of the receiving channels; 
b) Signal-to-Noise Ratio (SNR) of the event; c) Calculated phase difference between the 
received signals; d) Azimuth and elevation for each interpulse period in which the event was 
observed; e) Altitude of the detected meteor; f) Horizontal position of the detected meteor 


Figure 2: MIF modeled diurnal rate of meteoroid areal density at the MU radar geographical 
location during the observation times for the meteoroid mass range 0.1-100 pg 


Figure 3: Expected radiant distribution observable from the MU radar location for each of the 
experiment dates and meteoroid mass range 0. 1-100 pg, as modeled by the MIF 


Figure 4: Expected incoming speed distribution over the MU radar location for each of the 
experiment dates and meteoroid mass range 0. 1-100 pg, as modeled by the MIF 


Figure 5: Speed v. Altitude distribution of the MU radar (46.5 MHz) meteor observations 



Figure 6: RCS v. Altitude distribution of the MU radar (46.5 MHz) meteor observations 


Figure 7: Meteor RCS v. ablation altitude for various meteoroid masses and speeds, modeled 
for several HPLA radar systems with the RCS sensitivity threshold shown 


Figure 8: Observed distribution of the angle between the meteoroid velocity vector and local 
vertical 


Figure 9: Modeled and observed rates of meteor areal density corresponding to the MU 
observation periods 


Figure 10: Modeled and observed meteor speed distributions derived from the MU meteor 
observations 


Figure 1 1 : Modeled meteor radiant distributions predicted to be observed by the MU radar 


Figure 12: Observed meteor radiant distributions derived from the MU meteor observations 



